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Investigation of the non-uniform flow angularity effects on the Ares I DAC-1 in the 
Langley Unitary Plan Wind Tunnel are explored through simulations by OVERFLOW. 
Verification of the wind tunnel results are needed to ensure that the standard wind tunnel 
calibration procedures for large models are valid. The expectation is that the systematic 
error can be quantified, and thus be used to correct the wind tunnel data. The corrected 
wind tunnel data can thenn be used to quantify the CFD uncertainties. 


Nomenclature 

a angle of attack, degrees (°) 

a up flow angularity correction on angle of attack, degrees (°) 

AOA angle of attack 

Cn force coefficient in normal direction 

Ca force coefficient in axial direction 

Cmy pitching moment coefficient in y direction 

C t un coefficients corresponding to tunnel runs 

C n t coefficients corresponding to free air runs 

A re f reference area 

CFD computational fluid dynamics 

Cn normal force coefficient (Fn / qA re f) 

OVERFLOW Navier-Stokes computational fluid dynamics solver 
UPWT Unitary Plan Wind Tunnel (Langley) 


Note: To make this document ITAR free, ITAR sensitive figures have had their scales removed, and data 
tables have had their values scaled by an arbitrary number. 

I. Introduction 

T he Langley Unitary Plan Wind Tunnel, when operating at higher Mach numbers, exhibits an flow 
angularity angle upwards of 1.5°, according to the tunnel calibration report. * 1 As a result, large models 
such as the 0.01-scale model of the Ares I DAC1 experience a considerable amount of (potentially non- 
uniform) flow angularity over the length of the model. This can result in a different angle of attack (AOA) 
offset than determined by the standard wind tunnel data reduction procedures, introducing an unaccounted 
bias error into the data. Thus, verification of the wind tunnel results is needed to ensure that the standard 
calibration procedures are still appropriate for such large models and to quantify bias errors in the Ares 1 
DAC-1 data. 

OVERFLOW, a Navier-Stokes computational fluid dynamics (CFD) solver, can be used to validate 
the data obtained from the wind tunnel tests. Typically this is not done due to the unknown errors and 
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uncertainties in the CFD data. Traditionally wind tunnel data are used to validate and determine the 
uncertainty of the CFD data. The uncertainty and accuracy of the wind tunnel data are estimated by 
multiple wind tunnel runs and tunnel-to-tunnel comparisons. However since wind tunnel tests and numerical 
simulations of the full Navier-Stokes equations are expensive, an attempt should be made to extract as much 
information as possible and to assess the accuracy and uncertainty of both. In this study CFD is used to 
assess the accuracy and uncertainty of the wind tunnel data by avoiding the issue of unknown accuracy and 
uncertainty of the CFD data. The expectation is that much of the systematic error can be quantified and 
thus can be used to correct the wind tunnel data. The corrected wind tunnel data can then be used to 
quantify the CFD error and uncertainties. 

In this study, simulations of a simplified model (Ares I DAC-1) in the presence of simulated wind tunnel 
flow are performed at select angles of attack to determine the flow angularity effects on the AOA and other 
flow non-uniformities in the test section, to examine how well the solution performs after correcting for the 
effective flow angularities, and to verify the accuracy of the correction. 

II. Overview 

The goal of the present work is to run simulations of the Ares I DAC-1 model at 4 different angles of 
attack, simulating the wind tunnel tests. Simulations of the model in the test section of the wind tunnel 
are run at ±2° to determine the effective flow angularities by finding the AOA where Cjy= 0. The model 
is also run at ±7°+a up . These are called the correct wind tunnel simulations. The corrected results will be 
compared to the +7° free air run. The free air cases are considered to be the ’’true” results and the deviation 
of the corrected wind tunnel simulations from the free air simulations are due to the flow angularity and other 
test section flow non-uniformities. The data reduction procedure used for the wind tunnel data is applied 
to the wind tunnel simulation data. The values of the deviations will show the accuracy of the corrections. 
All of the numerical results reported in this work are produced with OVERFLOW 2.0aa. 2,3 

The entry section of the UPWT is simulated to obtain a much more representative flow state at the 
inlet station of the test section. Inflow conditions into the throat are set such that the static/stagnation 
conditions match the calibration report. For simulations involving test section 2, some smoothing is applied 
to the upper wall to aid convergence. 


III. Geometry 

The model tested is the 0.01-scale DAC1 as shown in Figure 1. However, all the values, including the 
Reynolds number, for the OVERFLOW simulation are scaled up to full scale. The reason for doing this is 
the original computations for the DAC1 configuration are performed at full scale, and that same grid system 
is used for the current simulations. The vehicle is placed such that the tip of the abort tower is at x=916.973 
in. when at 0° AOA, which translates to about 0.76 ft downstream of the test section inlet (x=0.0 in.) in the 
UPWT. The geometry (including the rocket and sting) is inserted into the tunnel, and the point of rotation 
is held constant ( ^rotation— 5348.186 in.). The test section modeled in the simulation is also scaled to 1, 
making the 4x4 test section 400x400 ft. Figure 2 shows the numerical model of the vehicle inside the tunnel, 
and Figure 3 shows a longitudinal (in the symmetry plane) cut. 

The model used for this simulation is the DAC1 vehicle analyzed previously. It is an axisymmetric body 
with the stiffener and kick rings on the booster as well as the LAS motors on the escape tower removed. The 
grid system consists of five zones. The first is a cap grid starting from the very tip of the rocket extending 
upstream with a very high resolution to resolve the stagnation point at the tip. The second grid is the tower 
grid system. The third grid includes the overall body of the rocket up to the end of the first solid rocket 
booster. The fourth grid includes the tail end of the booster with the aft skirt. The fifth grid is a simple 
cylinder modeling the sting. All of these grids are stretched grids with tight spacing near the body and 
coarsening as it extents radially outward. The test section grids are modeled using data from Reference 1, 
and these grids consist of two zones, a wall-normal stretched polar grid extending from the walls into the 
test section, and a core grid to remove the polar singularity of the wall grid. Table 1 outlines the number of 
grid points in each zone. 
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Figure 1: DAC1 Configuration shown with the model coordinate system. The tunnel coordinate system 
origin is 916 in. upstream of the model coordinate system origin 



Figure 2: Tunnel Test Section & Rocket Configuration 
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Figure 3: Tunnel Test Section & Test Article Symmetry Plane View 


Table 1: Number of Grid Points Per Zone 


Zone 

jt Points 

Cap 

80053 

Tower 

830428 

Body 

6113637 

Skirt 

2140506 

Wall 

9090000 

Core 

1440000 

Total 

20764877 


IV. OVERFLOW Simulation 

OVERFLOW 2.0aa is used for the CFD simulations. 2,3 The Spalart-Allmaras turbulence model is used 
to simulate the turbulence. 4 The distance function used in the turbulence model is modified to exclude 
the minimum wall distance, because it causes an error in the calculation of the eddy viscosity, causing the 
boundary layer to grow too fast resulting in an oblique shock forming off the inlet walls. Full multi-grid cycles 
are used to accelerate the initial system convergence, and additional iterations are performed afterward to 
further converge the solution. 

V. Calculating Flow Angularity Effects 

The flow angularities caused by the growing boundary layer and the changing area of the tunnel entry 
section will produce an effective change in the AO A for the overall flow field. To calculate this effective flow 
angularity, or a up , a procedure similar to the wind tunnel data reduction reported in Reference 5, is used. 
Runs are performed at ±2°, and the coefficients of normal force are compared. Using the AOA and Cn from 
the two runs, a line is fitted to these points and a up (i.e. where CA=0) is calculated. Equation 1 shows 
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how a up is calculated, and Table 2 shows the results of the ±2° simulation and the calculated a up . Mach 
and flow angularity contour plots of the runs performed and their corresponding residuals can be found in 
Appendix A. Note that a negative a up indicates a positive AOA for the flow angularity. 


OL U p = «2° - Cn 2 o («2° - ol-2°) /(Cn 2 o ~ Cn _ 2 o) (1) 

Table 2: Calculating Flow Angularity Correction 


AOA 

Cn 

2 

0.2698 

-2 

-0.2069 

&up 

-0.2639 


VI. Calculating the Correction for the Pitch Moment 

Due to model and flow field asymmetries the pitch moment is not necessarily zero when Cn= 0. The 
Cn and Ca are corrected by the angle of attack adjustment with the flow angularity. The same flow angularity 
will not always correct the pitch moment. Therefore the pitch moment is corrected by shifting the entire 
C my ~ a curve by A C my such that C my = 0 at a = a up . The shift is determined by 


A C m y,a up — Cmy ? 2° — {C my , 2° ~ C m y,-2° ) {OL up ~ a_ 2°)/( a 2° ~ «-2° ) (2) 

The definitions of the flow angularity and the C my shift are shown graphically in Figure 4. Again note 
that a negative value for a up indicates a positive AOA for the flow angularity. 


5 of 31 


American Institute of Aeronautics and Astronautics 



Figure 4: Correction Factors for Flow Angularity and AC my Shift 


VII. Correcting for the discrepancies between the inlet Mach number and 

Pressure and the nominal values 

Since the Mach numbers and pressures at the entrance of the test section are not quite equal to the nominal 
values, the force and moment coefficients have been rescaled. The same procedure used in Reference 1 to 
obtain the calibrated Mach number and pressure are followed. To determine the calibrated Mach number 
and pressure, the values at select locations in the x-z plane of the test section are averaged. Five longitudinal 
(x-dir) survey stations are used between 0 in and 48 in, and four vertical (z-dir) survey stations are used at 
±10.8in and at ±3.6in. For the vertical direction, the two probes closest to the center-line are given full 
weight, while the outer two data points given half weights. For test section 1, the first and last longitudinal 
points are given half weights and the rest given full weights. For test section 2, all longitudinal points are 
weighted equally. 

Because of the difference between the nominal values (free stream) and the calibrated values, the F&M 
results have to be rescaled by (P n om/ Piniet) * (M nom /Mi n i e t ) 2 . For example, for test section 1 at nominal 
Mach number of 1.6, the calibrated Mach number is 1.588. The corresponding nominal and inlet pressures 
are 1/1.4 and 0.71921, respectively. 
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VIII. Evaluating the correction for the effective flow angularity 

To evaluate how well the calculated a up predicts the effective flow angularity, runs at the corrected angles 
for ±7° are performed, and the normal forces are analyzed. The corrected ±7° AOA is defined as ±7° + 
a up . The same initial conditions applied to the ±2° runs are used. Additionally, these runs are compared 
to a free air run at +7°. Contour plots of the Mach numbers and flow angularities and their corresponding 
residuals can be found in Appendix A. 

IX. Calculating Boundary Layer 

The boundary layer thickness, displacement thickness, and momentum thickness of the tunnel simulations 
are also calculated. The boundary-layer thickness is defined as the distance from the wall to the point where 
u(y) = 0.99 u 0 , where u Q is taken to be the streamwise component of velocity measured at the centerline 
of the entrance to the test section (at 0,0,0). The compressible definition for displacement and momentum 
thickness is used. Displacement thickness is taken to be 



and momentum thickness is taken to be 



where p Q and u Q are the values at the centerline of the entrance to the test section (at 0,0,0). Both equations 
are integrated numerically using the trapezoidal rule. 

X. Results of Simulating Complete UPWT 

The difference in the results shows that the correction applied does not account for all of the flow 
angularities or Mach non-uniformity effects. It is very likely that interactions between the tunnel walls and 
test section inflow conditions and model are not adequately compensated for by the correction. Because the 
model is quite large in comparison to the tunnel test section cross-sectional area, tunnel-model interactions 
are more likely to occur. 

Figure 5 shows the complete tunnel geometry for test section 1. Figure 6 shows the grid system and 
Figure 7 shows the numerical solution of the empty tunnel in terms of Mach contours. Figures 8 and 9 show 
the Mach contours in the symmetry plane of the test section obtained from the numerical simulations and 
the experimental data. Unfortunately the experimental data are not taken at exactly the desired inlet Mach 
number. The closest Mach number to the simulation Mach number of 1.6 obtainable from the calibration 
report is Mach 1.469. Figures 10 and 11 show the corresponding flow angularities from the simulations and 
experiment. 

Similar simulations are performed at Mach = 2.5 for the low speed test section and at Mach 2.5 and 4.0 
for test section 2. Figure 12 shows the complete tunnel geometry for test section 2, and figure 13 shows the 
grid system. A similar sequence of results as shown above are presented in Figures 14 through 18 for the 
Mach 4 case. 

The boundary layer thicknesses at the entrance of the test section are also in good agreement with the 
tunnel data (Figure 19). The trends of the data are similar to the tunnel data, while the overall shift in the 
data could be a result of the simulation model. A thinner boundary layer is expected, because only the inlet 
through the test section are included in this simulation. 

The same procedure as described above is followed to compute all 4 cases and the discrepancies between 
the computed tunnel results at the corrected angles of attack and pitch moment shifts and the free air (no 
tunnel) simulations are determined. These results are shown in Figure 20 and 21. Note: Figure 20’s values 
for normal, axial, and pitch coefficient are rescaled by arbitrary numbers to make it ITAR free. 

With the more representative entry section the inlet flow conditions of the test section obtained from 
the simulations are in much better agreement with the reported results of the calibration report. 1 While the 
Mach contours and the flow angularities are not in perfect agreement with the calibration report, they are 

7 of 31 


American Institute of Aeronautics and Astronautics 


representative and should provide a good first order indication of the errors in the reported wind tunnel data 
due to flow non-uniformities (Mach and flow angularities). 
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Figure 5: UPWT Entry and Test Section 1 



Figure 6: Grid System for Entire UPWT with Test Section 1 at Mach 1.6 



Figure 7: Computed Mach Contours of Entire UPWT with Test Section 1 at Mach 1.6 
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Figure 8: Computed Mach Contours of UPWT Test Section 1 at Mach 1.6 
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Figure 9: Experimental Mach Contours of UPWT Test Section 1 at Mach 1.469 
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Figure 10: Computed Flow Angularity of UPWT Test Section 1 at Mach 1.6 


in. y=0, Oin. (0. 0cm) 



_i i i i 1 — 1 1 ' 1 . 

40 60 SO 100 120 140 160 ISO 200 cm 

Longitudinal distance, x 


Figure 11: Experimental Flow Angularity of UPWT Test Section 1 at Mach 1.57 
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Figure 12: UPWT Entry and Test Section 2 



Figure 13: Grid System for Entire UPWT with Test Section 2 at Mach 4.0 
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Figure 14: Computed Mach Contours of Entire UPWT with Test Section 2 at Mach 4.0 
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Figure 15: Computed Mach Contours of UPWT Test Section 2 at Mach 4.0 
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Figure 16: Experimental Mach Contours of UPWT Test Section 2 at Mach 3.953 
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Figure 17: Computed Flow Angularity of UPWT Test Section 2 at Mach 4.0 
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Figure 18: Experimental Flow Angularity of UPWT Test Section 2 at Mach 4.185 
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Figure 19: 


Boundary Layer Properties of UPWT and Overflow Computational Results 
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Figure 20: Wind Tunnel Non-uniformity Induced Errors in the Longitudinal F and M Coefficients for 

DAC1 in the UPWT Test Sections 1 and 2 
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Induced Errors 



Figure 21: Wind Tunnel Non-uniformity Induced Errors in the Longitudinal F and M Coefficients for 

DAC1 in the UPWT Test Sections 1 and 2 


As shown by Figure 20 and 21, the errors of F&M data are still substantial. The error is defined as 
( C tU n ~ C nt )/C nt , where C is any one the force or moment coefficients. Figure 21 shows the percent error 
for the nominal +7° angle of attack. For Cn the errors range from 1.5% at Mach 1.6 to 4% at Mach 4. The 
Ca errors are 0.2% at Mach 1.6 to ~ 0.4% at Mach 4. The corresponding pitch moment errors range from 
~ 0% to ~ 4.5% over the same Mach number range. 

These errors can now be used to correct the reported wind tunnel data. 5,6 The rationale behind this is 
quite simple. The same grid system around the DAC1 configuration, the same turbulence model, flow solver, 
boundary conditions, and convergence levels are used for all five cases. Therefore the effects of the CFD 
uncertainties are constant for all five cases and effectively vanish when increments are computed from these 
five cases. These increments (errors) can be used to correct the reported wind tunnel data without any fear 
of contamination from the CFD uncertainties. Any other effects in the CFD data will be due to the flow 
non-linearities occurring between a and a = a up . 

Following the guidelines of Reference 5, the corrected values are not used for the Mach 2.5 in Test 
Section 1. The reason for excluding this case is also evident from Figure 20. For this case the effective flow 
angularity is -0.960. The magnitude of both is substantially larger than for the other 3 cases. Using only the 
three remaining cases the corrected values are shown in the following figures along with the original wind 
tunnel data and the prior CFD results reported at the 2nd CLV Aero TIM meeting. The corrections are 
only developed for Mach 1.6, 2.5, and 4.0 at a = +7°. The axial force, normal force, and pitch moment 
coefficients are shown in Figures 22, 23, and 24, respectively. The corrected wind tunnel data is in closer 
agreement with the CFD results. The errors in Cn and C my are now less than 2% at Mach 2.5 and 4 and less 
than 5% at Mach 1.6. The discrepancies between the CFD and the uncorrected wind tunnel are in excess 
of 5%. The errors are still fairly large at Mach 1.6. This may be due to the larger tunnel wall interference 
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effects at the lower Mach numbers. These wall effects are not likely to occur at the higher Mach numbers. 
It is also possible that the test section inlet flow angularities and flow non-uniformities are not resolved with 
enough precision. These discrepancies need further scrutiny. 
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Figure 22: Comparison of the Axial Force Coefficients for DAC1 in the UPWT Test Sections 1 and 2 with 

the corrected UPWT data and the Overflow results 
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Figure 23: Comparison of the Normal Force Coefficients for DAC-1 in the UPWT Test Sections 1 and 2 
with the corrected UPWT data and the Overflow results 
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Figure 24: Comparison of the Pitch Moment Coefficients for DAC-1 in the UPWT Test Sections 1 and 2 
with the corrected UPWT data and the Overflow results 


XI. Grid Refinement 

To see if refining the grid would yield a different solution, simulations are run using an adjoint-based 
refined grid (from Cart3D). 7 Figure 25 shows the Cart3D adjoint-based refined grid, as well as the corre- 
sponding domain of influence and the domain of dependence regions. The near body grids are refined using 
the same near-body spacing as the adjoint-based refined grid. Additionally, an off-body grid is added to 
increase the grid resolution encapsulated by the domain of influence and domain of dependence regions. The 
runs with the new adjoint grid are performed at Mach 1.6, and compared with the old grid system. Figure 26 
shows the adjoint-based grid that is added. Table 3 outlines the number of grid points per zone. Figure 29 
shows the induced errors for the adjoint grid, compared with the original grid. The difference in the induced 
errors between the two different grid systems is small, implying that the adjoint grid did little to change the 
solution. Note: Figure 29’s values for normal, axial, and pitch coefficient are rescaled by arbitrary numbers 
to make it ITAR free. 



Figure 25: Adjoint-based grid from Cart3D 
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Figure 26: Grid for the adjoint-based body grids for the DAC-1 model 


Table 3: Number of Grid Points Per Zone, Adjoint Grid 


Zone 

jt Points 

Cap 

80053 

Tower 

1447638 

Body 

6817908 

Skirt 

3408954 

Near-Body 

4538575 

Wall 

9090000 

Core 

1440000 

Total 

46296194 


XII. Comparison with Baldwin-Barth turbulence model 

To determine the effect of the turbulence model on the results of the simulations, runs are performed 
again at M = 1.6 using the Baldwin-Barth turbulence model. 8 Figures 27 and 28 show the Mach and 
flow angularity contours in the symmetry plane of the test section obtained from using the Baldwin-Barth 
turbulence model. Again the same procedure is used to determine the effective flow angularity, and evaluate 
how well a up performed. The adjoint-based grids are used in this simulation. Additionally, for consistency, 
the free-air case is run with the same turbulence model (Baldwin-Barth). The same normalizations are 
applied to the results, which are shown in figure 29. It seems with the thinner boundary layer (figure 19, the 
non-uniform flow angularity effects has a larger impact on axial force and pitching moment, as the errors 
are 2.5x greater than the SA turbulence model results. Note: figure 29’s values for normal, axial, and pitch 
coefficient are rescaled by arbitrary numbers to make it ITAR free. 
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Figure 27: Mach Contours of UPWT Test Section 1 at Mach 1.6, Baldwin-Barth Turbulence Model 



Figure 28: Computed Flow Angularity of UPWT Test Section 1 at Mach 1.6, Baldwin-Barth Turbulence 

Model 
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Figure 29: Comparison of Induced Errors for Mach 1.6 configuration, original grid, adjoint grid, and 

different turbulence models. 


XIII. Summary of Procedure 

This section gives a brief outline of the procedure used to correct for the tunnel effects on the reported 
wind tunnel data. An attempt is made to follow the standard wind tunnel data reduction technique used to 
eliminate test section inlet flow non-uniformities and flow angularity effects on the angle of attack. 

1. Compute grid system in free air at a = a nom for clean DAC1 configuration. 

2. Compute solution in free air at a = a nom . 

3. Compute grid system for empty tunnel with entire inlet section and test section. 

4. Compute solution in empty tunnel. 

5. Compare flow non-uniformities and boundary layer thickness with experimental data from Calibration 
report 1 . 
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6. Install free air grid system in tunnel at a = ±2°. 

7. Determine effective (or integrated) flow angularity (a up ). 

8. Install free air grid system in tunnel at a ± a up . 

9. Compare integrated forces ( Cn ) at a = ±7 ± a up with free air forces at a = 7°. 

10. Differences are the wind tunnel errors in Cn due to the test section flow non-uniformities and tunnel 
wall effects. 

11. Correct wind tunnel data in CV for the wind tunnel errors. 

12. C my may not be zero at a = a up (flow angularity correction for Cn)- Compute the additional correction 
to Cmy to force (C my ± AC my ) = 0 at a = a up . 

13. Correct wind tunnel data in C my (ie. shift by AC my ). 

14. No additional correction beyond the a up is applied to Ca- 

15. Do not correct for side forces (C y ), no way of doing this. 

16. Use corrected wind tunnel data to validate the CFD solutions. 

17. Rationale: 

(a) Same grid system around the DAC1 configuration and same turbulence model, flow solver, bound- 
ary conditions, and convergence levels are used for all five runs for each case. 

(b) Because of (a), the effects of the CFD uncertainties are removed and the increments (errors) can 
be used to correct the reported wind tunnel data. Any other effects in the CFD data will be due 
to the flow non-linearities between a and a ± a up . 

18. Add uncertainty bounds due to the tunnel Mach number uncertainties. For example the calibration 
report indicates that the uncertainty of the measured Mach number is ±0.01. 

XIV. Summary 

This procedure has shown that it is possible to remove some of the systematic uncertainty from the wind 
tunnel data by the use of numerical simulations to compute the discrepancies between the tunnel results and 
free air results. The technique is fairly compute-intensive, requiring four complete simulations for each set 
of flow parameters. However only a limited number of cases need to be corrected to show the extent of the 
required corrections. In the example shown above the systematic error or discrepancy between the CFD and 
wind tunnel data has been reduced from more than 5% to less than 3%. There is still an error larger than 
3% at Mach 1.6. This Mach number needs further investigation. The accuracy of using numerical simulation 
to capture the test section inflow state needs to be ascertained. Unfortunately the calibration report has 
limited data on the inlet flow state. 1 

With more of the tunnel artifacts removed from the wind tunnel data, the corrected wind tunnel data 
are a much more appropriate experimental data set for validating CFD data. Typically the CFD data is for 
free- air flow fields and wind tunnel data are not. Much of the systematic errors in wind tunnel data are due 
tunnel effects, either walls or non-uniform inlet flow fields. Current wind tunnel data reduction techniques 
do not seem to eliminate all of the tunnel artifacts from the data. The proposed technique provides a method 
to further reduce tunnel artifacts from the wind tunnel data. 

The method needs to be tested for other force and moment coefficients, in particular the roll moment 
coefficient for CLV configurations with protuberances. The roll moments are believed to be especially 
sensitive to tunnel artifacts such as flow non-uniformities and variable flow angularity angles. 
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Appendix 



Figure 30: Mach Contours for Tunnel 1 Mach 1.6 configuration at +7° AO A 



Figure 31: Flow Angularity Contours for Tunnel 1 Mach 1.6 configuration at +7° AO A 
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Figure 32: Mach Contours for Mach 1.6 Free Air configuration at +7° AO A 



Figure 33: Flow Angularity Contours for Mach 1.6 Free Air configuration at +7° AO A 
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Figure 34: Mach Contours for Tunnel 2 Mach 4 configuration at +7° AOA 



Figure 35: Flow Angularity Contours for Tunnel 2 Mach 4 configuration at +7° AOA 
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Figure 36: Mach Contours for Mach 4.0 Free Air configuration at +7° AOA 



Figure 37: Flow Angularity Contours for Mach 4.0 Free Air configuration at +7° AOA 
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Figure 38: Mach Contours for Tunnel 1 Mach 1.6 configuration at +7° AOA, Adjoint-Based Grid 



mill 


Figure 39: Flow Angularity Contours for Tunnel 1 Mach 1.6 configuration at +7° AOA, Adjoint-Based Grid 
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Figure 40: Mach Contours for Tunnel 1 Mach 1.6 configuration at +7° AOA, Baldwin-Barth Turbulence 

Model 



Figure 41: Flow Angularity Contours for Tunnel 1 Mach 1.6 configuration at +7° AOA, Baldwin-Barth 

Turbulence Model 
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Mach Number: 1.5 1.59 1.68 



Figure 42: Mach Contours for Mach 1.6 Free Air configuration at +7° AOA, Baldwin- Barth Turbulence 

Model 



Figure 43: Flow Angularity Contours for Mach 1.6 Free Air configuration at +7° AOA, Baldwin-Barth 

Turbulence Model 
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Residual History 



Time Step Number 


Figure 44: Typical Residual for Full Tunnel simulation 


Residual History 

LoglO(L2 norm of RHS) 



Figure 45: Typical Residual for +7 AO A Tunnel simulation 
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